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ABSTRACT 

We construct a sample of extremely red objects (EROs) within the UKIDSS Ultra Deep Sur- 
■ vey by combining the Early Data Release with optical data from the Snbam/XMM-Newton 

Qh Deep Field. We find a total of 3715 objects over 2013arcmin 2 with R — K > 5.3 and 

K < 20.3, which is a higher surface density than found by previous studies. This is partly 
£^ ■ due to our ability to use a small aperture in which to measure colours, but is also the result 

^| of a genuine overdensity of objects compared to other fields. We separate our sample into 

C$ passively-evolving and dusty star-forming galaxies using their RJK colours and investigate 

their radio properties using a deep radio map. The dusty population has a higher fraction of 
individually-detected radio sources and a higher mean radio flux density among the undetected 
. objects, but the passive population has a higher fraction of bright radio sources, suggesting 

' that AGNs are more prevalent among the passive ERO population. 
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1 INTRODUCTION 

Extremely Red Objects (EROs) are variously defined on the basis 
of an optical-near-infrared colour, usually R — K or / — K. These 
red colours can be due either to an evolved, passive stellar popula- 
tion at z > 1, or the presence of large quantities of dust associated 
with starburst activity (although edge-on discs may also contribute 
significantly; Yan & Thompson 2003). While the former population 
reveals the stellar mass already assembled in galaxies as a result of 
star formation at early cosmic times, the latter population reveals 
the star formation at z ~ 1 which, due to its dust-enshrouded na- 
ture, may be missed by other surveys. Since these two populations 
represent very different routes to the creation of present-day galax- 
ies, there is a clear desire to reliably separate them and use them as 
tests of galaxy formation models. 

Existing samples of EROs (for which we adopt the most com- 
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mon definition of R — K > 5.3) have been relatively small (typ- 
ically no more than about 100 objects), due to the difficulty in ob- 
taining near-infrared imaging data with the necessary depth and 
areal coverage (e.g., Daddi et al. 2000; Smail et al. 2002; Mous- 
takas et al. 2004; Georgakakis et al. 2006). The arrival of infrared 
cameras with large fields of view, such as WFCAM (Casali et al., 
in preparation), has largely eliminated this problem, allowing sam- 
ples of EROs to be constructed which are large enough to permit 
further separation by magnitude, colour, and other properties. Fur- 
thermore, by locating infrared surveys in well-observed patches of 
sky, it becomes possible to study the multi-wavelength properties 
of EROs and learn more about their nature. 

In this Letter, we combine near-infrared data from the 
UKIDSS Ultra Deep Survey Early Data Release (Dye et al. 2006; 
Lawrence et al. 2006) with optical data from the Subaru/XMM- 
Newton Deep Field (SXDF; Furusawa et al., in preparation) to con- 
struct a sample of EROs to faint magnitudes {K > 20) over an 
area larger than half a square degree. The Ultra Deep Survey will 
eventually cover 0.8 deg 2 to 5a point source limits of J = 25.0, 
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H = 24.0, and K = 23.0 and has been sited in the SXDF be- 
cause of the excellent deep multi-wavelength data in this region. 
We determine photometric redshifts for all objects and separate 
them into dusty star-forming and passively-evolving galaxies on 
the basis of their RJK colours. We use deep radio data to deter- 
mine the typical radio flux densities and star formation rates as a 
function of colour and magnitude. A more detailed analysis, in- 
cluding clustering and X-ray properties, will be undertaken after 
the UKIDSS DR1 release, which is expected to be a magnitude 
deeper than the present data, and will provide more reliable pho- 
tometry. All magnitudes in this paper are on the Vega system. The 
SXDF catalogue is calibrated on the AB system, and we use the 
conversion i?v cg a = Rab — 0.219 (H. Furusawa, private commu- 
nication). We adopt a cosmology with Hq — 70kms _1 Mpc -1 
and fi m = 1 — S1a = 0.3. 



2 THE ERO SAMPLE 

We start with the improved version of the UDS EDR if-band cat- 
alogue produced by Foucaud et al. (2006), and determine the com- 
pleteness for galaxies as a function of if -band magnitude. Galaxies 
with K ps 20 have a median FWHM = 1.2" and so we added ob- 
jects of this size to the if -band WFCAM images and measured the 
fraction that were recovered by the extraction process. The 50% 
completeness limit varies slightly between detectors, but averages 
to K = 20.3. We therefore eliminated all objects from the cata- 
logue of Foucaud et al. (2006) with (total) magnitudes fainter than 
this limit before proceeding further. Our magnitude limit is brighter 
than that of Foucaud et al. since they calculate the completeness for 
unresolved sources. The optical data reach R ~< 26.7 (5a, <f>2"). 



2.1 Optical-infrared matching 

In order to facilitate matching of the near-infrared catalogue with 
objects in the optical imaging, objects in various regions of sky 
were next removed from the K < 20.3 catalogue. Only the cen- 
tral region of the UDS image, with uniform noise, was included in 
the matching process. Objects were eliminated whose UDS coor- 
dinates placed them outside the optical imaging region, or within 
2 arcsec of the edge, to ensure that the photometric aperture of any 
optical match was entirely within the image. Finally, objects which 
lay within the halo or CCD bleed of a bright star were removed 
from the catalogue. The remaining effective area is 2013 arcmin 2 . 

The UDS catalogue was first matched to the catalogues of Fu- 
rusawa et al. (in preparation), with the closest match within 1 arcsec 
being used. However, visual inspection of the unmatched objects 
revealed that many pairs of objects remain blended in the opti- 
cal catalogues, so the UDS catalogue was instead matched to the 
catalogues produced by Simpson et al. (2006) with a more severe 
deblending threshold. The i'-band catalogue was used first (as the 
i-band data are the deepest), followed by the z', R, and B-band cat- 
alogues. This reduced the number of unmatched objects by about 
half, to 426. These were all inspected by eye, and the vast major- 
ity (301) were discovered to be artifacts such as cross-talk, satel- 
lite trails, or array defects not removed in the reduction process 
(see Dye et al. 2006). These were deleted from the catalogue and 
not considered further. Most of the remaining 125 sources are still 
blended in the optical catalogue, although some lack any optical 
counterpart. Optical photometry was determined for these objects 
directly from the Suprime-Cam images using the UDS catalogue 
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Figure 1. Plot of r_2 (Kron 1980) against J — K colour for all objects 
with K < 20.3. The solid line shows the criterion used to separate stars 
from galaxies. 

coordinates. A total of 443 objects have measured i?-band aperture 
magnitudes fainter than the 5a limit of the optical data. 

In what follows, colours are measured in a 2-arcsecond aper- 
ture, while magnitudes are total magnitudes from SExtractor's 
MAG_AUTO parameter (Bertin & Amouts 1996). The if -band im- 
age has a narrower psf than the other images (FWHM=0.69 arcsec 
compared to 0.81 arcsec), and we address this point below. 

2.2 Star-galaxy separation 

We classify objects as stars or galaxies by the same method as 
Simpson & Rawlings (2002), using a combination of J — if 
colour and Kron (1980) radius, r_2. There is a clear bimodality 
in this plot (Fig.0 because galaxies rapidly become redder than 
stars with increasing redshift, while those with similar colours to 
stars are nearby and hence clearly extended. While there was a 
clear separation between the two groups in the small sample of 
Simpson & Rawlings (2002), the much larger number of objects 
in the UDS EDR means that many are scattered into 'no man's 
land', and we determine the locus of the minimum surface den- 
sity of objects in Fig.Qas our separation criterion. All objects with 
r_2 < 0.41 — 0.1( J — if) are classified as stars. However, the com- 
bination of depth and area in the UDS will result in the detection of 
a number of very late-type stars with extremely red J — if colours, 
and these will not be excluded by this single criterion. We therefore 
also classify as stars any objects with r_2 < 0.3 and K < 19.5 
(morphological classification is not possible fainter than this). 

Our galaxy sample will suffer contamination from L dwarf 
stars, which have J — K > 1 (e.g., Hewett et al. 2006). Objects 
with J — K > 1.3 and if > 19.5 will be classified as galaxies, 
and these red colours are possessed by stars of spectral types Ll- 
L8 (Table 9 of Knapp et al. 2004), whose surface density to if < 
14.5 is 0.038 deg" 2 (Reid et al. 1999). If these are isotropically 
distributed on the sky (so the number expected is proportional to 
the Euclidian volume sampled), we expect ~ 60 objects with K < 
20.3 in our survey area. These are therefore a small contaminant. 

2.3 Number counts 

To construct our ERO sample, we first investigate the effect of 
the superior if -band image quality by Gaussian-smoothing objects 
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Figure 2. Differential source counts for objects in the SXDF/UDS com- 
pared to a selection of results from the literature. 



with R — K > 5.3 to the same FWHM as the other images and 
measuring the flux loss in a 2-arcsec aperture. We find a median 
reduction in brightness of 0.04 mag, and apply this reduction to all 
A'-band magnitudes in the UDS EDR catalogue before construct- 
ing our ERO sample. In total we find 3715 (5128; 1464) objects 
with R - K > 5.3 (R - K > 5; R - K > 6). 

In Fig.|2|we compare our galaxy and ERO number counts with 
determinations from the literature. Our total galaxy number counts 
agree very well, except at the brightest magnitudes where there is a 
slight deficit due to the fact that we have excluded regions with very 
extended galaxies from our analysis. Our ERO counts are, however, 
~ 60 per cent higher than those of Daddi et al. (2000). This could 
be due in part to the different ii-band filters used in these stud- 
ies, since Daddi et al.'s filter has a long red tail with its 10 per cent 
point at 8200 A, while the Suprime-Cam filter has a sharper long- 
wavelength cutoff with a 10 per cent transmission point at 7200 A. 
Consequently, objects with a pronounced (4000-A) break between 
these two wavelengths will appear fainter in the Suprime-Cam im- 
ages, and hence redder in R — K. However, we note that we also 
find a much larger surface density of R — K > 5.3 objects than did 
Smail et al. (2002), who used the same instrument/filter combina- 
tion for their optical imaging. This may be partly due to the small 
aperture we use to determine colours, since contamination from 
nearby objects (which is obviously more likely for larger apertures) 
will almost certainly make the aperture colour bluer. We find that 
the number of EROs is decreased by 15-20 per cent (with a larger 
reduction for redder objects) if we use the colours measured in a 
4-arcsecond aperture. Finally, we note that our surface density is 
similar to that of Miyazaki et al. (2003) who studied a small re- 
gion within the SXDF, suggesting that this region is overdense at 
z ~ 1. For objects as clustered as EROs (ro ~ 12Mpc; Daddi et 
al. 2001), the effect of cosmic variance can be large, and we esti- 
mate a v ~ 0.2 using the prescription of Somerville et al. (2004). 
We stress that the depth of the 7?-band data (the 3<r limit in a 2- 
arcsecond aperture is R ~ 27 A; Furusawa et al., in preparation) 
means that photometric errors will not scatter a significant number 
of the more numerous bluer objects into the ERO class. 



3 NATURE OF THE ERO POPULATION 

Although the ERO population is composed of both passively- 
evolving red galaxies and dusty starburst galaxies, there does not 
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Figure 3. Histograms of the distribution of photometric redshifts, 2 p hot> 
for subsamples of EROs. In each plot, the solid lines show the distribution 
of bright (K < 19.5) EROs, while the dashed line show the distribution of 
faint (K > 19.5) objects. The numbers by the line-colour key indicate the 
total number of bright and faint objects in each class, respectively. 



appear to be a bimodality in any photometric property which al- 
lows these two classes to be readily separated. Pozzetti & Man- 
nucci (2000) used J-band photometry to attempt to separate the 
classes as heavily-reddened objects will not show the pronounced 
4000-A break apparent in the spectra of old galaxies, and hence 
will have bluer R — J and redder J — K colours. Smail et al. 
(2002) used a more sophisticated method of spectral template fit- 
ting with the code HyperZ (Bolzonella, Miralles, & Pello 2000) 
and classified objects according to whether the best-fitting model 
had a low or high extinction, finding that this produced broadly the 
same classification as the simpler method. Both methods rely on 
the accuracy of the photometry around A rcst ~ 4000 A, and our 
ability to classify EROs is therefore limited by the quality of the J- 
band photometry. Since the dividing line between the two classes 
of object is at J — K ~ 2 and the J-band photometry limit is 
J = 21.6 (5a, 4>2"), a classification is likely to be reliable only 
for the 1922 EROs with K < 19.5 and we do not believe a so- 
phisticated analysis is warranted by the present near-infrared data 
quality. We therefore opt to use Pozzetti & Mannucci's (2000) sim- 
ple RJK colour criterion to separate our sample into passive and 
dusty galaxies. We find that 62% of K < 19.5 EROs are clas- 
sified as passively-evolving, although this fraction increases when 
a redder colour cut is used, rising to 73% (77%) of objects with 
R — K > 6 (R — K > 6.5). These fractions are broadly consistent 
with, although perhaps slightly larger than, those determined from 
previous studies. 

We have computed photometric redshifts for all members of 
our ERO sample, in order that we can derive physical quantities 
from observed fluxes, by fitting synthetic galaxy templates to the 
BVRi'z'JK photometry, plus the two short- wavelength IRAC 
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bands (from the SWIRE survey of Lonsdale et al. 2003) where 
available, using a modified version of the HyperZ code (see Cira- 
suolo et al., in preparation, for a detailed description). We show 
the results of this analysis by plotting the distribution for vari- 
ous subsamples in Fig. [3] Objects whose best template fits have 
X 2 > 20.09 (corresponding to an unacceptable fit at 99 per cent 
confidence) are not plotted. 

Fig.|5|also reveals the predictable results that the reddest ob- 
jects are typically at higher redshift, as are the faintest objects. The 
redshift distribution for passive galaxies also appears to be slightly 
narrower and more strongly peaked than that for dusty ones. This is 
due to the passage of the 4000-A break through the filters, as galax- 
ies will only be red in R — K if the break is beyond the blue end of 
the R filter, while the Pozzetti & Mannucci (2000) criterion breaks 
down at z > 2 and classifies all EROs as dusty. We caution against 
a more detailed analysis of the different N(z) distributions given 
the significant difference between the photometric accuracy of the 
optical and near-infrared data, even for EROs. The deeper JK data 
from UKIDSS DR1 will provide better photometry, and hence a 
more reliable photometric analysis, for the present sample, as well 
as allowing a basic study of the fainter (K ~ 21) ERO population. 
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Figure 4. Radio properties of EROs with 18.5 < K < 19.5, R - K > 6, 
and classified as 'passive' . The histogram shows the radio image pixel val- 
ues at the locations of the sample objects (thick line), compared to the pixel 
values of the 2013 arcmin 2 over which EROs were selected (thin line). The 
inset shows the results of the image stacking, after rejecting individually 
detected sources. A clear positive signal can be seen in both analyses. 



4 RADIO PROPERTIES 

Simpson et al. (2006) have presented a 1 .4-GHz radio image of the 
SXDF taken with the NRAO Very Large Array which reaches an 
rms noise level of 12 (ijy beam -1 in the central regions. At the me- 
dian redshift of our sample, z ~ 1.3, a star formation rate (SFR) 
for M > 5Mq stars of IOMqvi- -1 produces a flux density of 

5 /iJy (Condon 1992). Although this is obviously below the detec- 
tion threshold for any individual source, we can make statistical 
detections of relatively modest SFRs in samples of > 100 objects 
since the 5" x 4" beam will encompass all radio emission from the 
host galaxy. We note that starburst galaxies will be unresolved at 
the resolution of this radio image. 

We study the mean radio properties of our sample in two ways. 
First, we produce small postage stamps from the radio map of 
Simpson et al. (2006) at our objects' positions and sum these af- 
ter weighting each by the inverse of the pixel noise variance at its 
position. We exclude all objects where any of the nine pixels (the 
scale is 1.25 arcsec pixel -1 ) in the radio image closest to the source 
position has a signal-to-noise ratio in excess of 3, so that the mean 
is not biased by a few objects associated with bright radio sources. 
We then fit an elliptical Gaussian profile to this stacked image us- 
ing Ivan Busko's n2gaussf it task in IRAF, and calculating the 
flux density from the fit using the prescription of Condon (1997). 
For the second method, we construct a histogram of pixel values 
from the radio map at the positions of objects in our sample, and 
fit a Gaussian to the pixels after they have been iteratively sigma- 
clipped to exclude values more than two standard deviations from 
the mean, thus automatically removing bright outliers. We refer to 
these two methods as the stacking and histogram methods, respec- 
tively. Since the histogram method fits a Gaussian to the convolu- 
tion of true source fluxes with the near-Gaussian pixel noise distri- 
bution, it is most reliable if the source flux density distribution is 
either itself Gaussian, or has a dispersion much less than that of the 
overall pixel noise. We prefer to use the image stacking method to 
derive mean radio flux densities, but confirm that the results from 
the histogram method are consistent in all cases. 

To determine the likelihood of spurious detections in our im- 
age stacks, we ran Monte-Carlo tests by producing 1000 lists of 



between 100 and 1000 sources each located at random positions 
within the unmasked region of the SXDF/UDS images. The mean 
flux density at the positions in each list was then calculated by both 
the stacking and histogram methods. Since the stacked radio im- 
ages are essentially signal-free, the Gaussian fitting technique was 
highly unstable and instead we measured the signal in a 5 x 5 pixel 
box in the centre of the image and corrected this for the 20 per 
cent flux loss expected for the beam size. Both methods produced 
a mean flux density of —0.2 /iJy, and we have therefore corrected 
our measured flux densities for this offset. The standard deviations 
on the mean for a sample of 100 sources were 0.9 and 2.1 /iJy for 
the stacking and histogram methods, respectively, and reduced with 
sample size N approximately as TV -1 / 2 . 

TableQshows the results of the stacking analysis for 18 ERO 
subsamples (separated by magnitude, colour, and classification), 
while Fig. [4] shows the stacking and histogram analyses for one 
subsample, indicating the reality of the statistical detection. From 
our analysis, galaxies classified as 'dusty' typically have twice the 
mean radio flux density of 'passive' galaxies, and the fraction of 
individually-detected objects is about twice as high (except in the 
faintest magnitude bin, where we have argued that the classifica- 
tion is unreliable). Smail et al. (2002) also find that most EROs 
with flux densities of a few tens of /iJy are dusty. However, while a 
smaller fraction of passive galaxies are detected individually, those 
that are tend to be brighter: in the 18.5 < K < 19.5 bin, 28 (13) 
of 104 individually-detected passive EROs have 5*1.4 > 100 ^tJy 
(Si. 4 > 300 /Jy), while only 10 (4) of 98 dusty EROs are as 
bright (for reference, the Fanaroff & Riley (1974) break is at 
5*1.4 ~ 6 mjy at z = 1.3). This suggests that a few per cent of the 
passive ERO population may be active galaxies (cf. Alexander et 
al. 2002), although such objects can also be found in the dusty pop- 
ulation (e.g., EDXS N2J21; Willott et al. 2003). Among the dusty 
population, we find that the SFR is uncorrected with magnitude, 
but may be correlated with colour. 

Our mean radio flux densities are consistent with those de- 
termined by Georgakakis et al. (2006), although our larger sample 
size produces much smaller uncertainties. We find that dusty galax- 
ies typically have about twice the 1 .4-GHz flux density of passive 
galaxies with the same colours and magnitudes, and hence have 
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Table 1. Average radio flux densities of ERO subsamples, determined from Gaussian fitting to the variance-weighted images (the 'stacking' method; see text). 
For each magnitude- and colour-selected subsample, the median photometric redshift is shown (only considering objects with x 2 < 20.09), together with 
the number of objects stacked and the number rejected from the stacking analysis due to being individually detected. The mean radio flux density and star 
formation rate (for stars with M > 5 Mq), adopting the median photometric redshift for the subsample, are also shown. 
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3.4 ±0.8 


12.2 ± 2.7 
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t Ellipse parameters fixed during fit due to high noise. 



about twice the star formation rate since the median redshifts of 
identically-selected passive and dusty EROs are similar (in comput- 
ing the median redshifts for our subsamples, we ignore cases where 
the best-fitting photometric redshift has \ 2 > 20.09). However, the 
true difference in star formation rates is likely to be larger than this 
since the individually detected dusty objects are likely to be the 
high end of a broad distribution of star formation rates whereas the 
individually detected passive objects are likely to be AGNs. As a 
final point, we note that the sample of 5834 non-ERO galaxies with 
18.5 < K < 19.5 and R - K < 5.3 has 440 objects detected at 
> 3cr significance, while the remaining sources have a mean flux 
density of 5*1.4 = 5.3 ± 0.2 /ijy, consistent with the value derived 
for the passive ERO population in this magnitude range. 



5 SUMMARY 

We have described the construction of a sample of extremely 
red objects (EROs) within the Subaiu/XMM-Newton Deep Field 
(SXDF) and UKIDSS Ultra Deep Survey (UDS) which is much 
larger than previous samples. We have measured a higher surface 
density than found in other fields, but comparable to that measured 
in a small sub-region of our field, suggesting the SXDF samples a 
cosmic overdensity at z ~ 1. After using RJK colours to separate 
the sample into passively-evolving (62% of EROs with K < 19.5) 
and dusty star-forming galaxies, we have used a deep radio map to 
study the radio properties of these two populations. Both popula- 
tions are detected at high significance in a stacked radio image, with 
the dusty sample having a mean flux density (and derived mean star 
formation rate) equal to approximately twice that of the passive 
sample. The mean radio flux density of the passive sample is con- 
sistent with that of non-ERO galaxies with the same magnitudes. 
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